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ABSTRACT 

o “f* o 4 1 

Interstellar lines arising in C , C , O , Si , 
and S + observed in the ultraviolet spectrum of C Oph 
by rocket spectrographic techniques are analyzed. The 
main results are that within a factor of 10 the 


abundances of C 


o *f* a 

O . and Si outside the H II region 


surrounding C Oph relative to the hydrogen abundance 

are equal to solar values. Hie lines in C° and S + imply 

that the interstellar matter is distributed among 

several clouds as indicated by high resolution visible 

spectra. The 2 P°^ 2 (°° 0079 eV ) fine structure level of 

4. 

the ground state configuration in C is strongly 
populated, the abundance ratio of the 2 j> 3/2 * eve * _to the 
2p i/2 S round level being 1.5. It is suggested that 
the excited C + ions are inside the Stromgren sphere 

-3 

where proton densities equal to or greater than 220 cm 
can collisionally excite the ions at sufficient rates. 
Stellar absorption lines of C TV (1548.2, 1550. 8$) and 
N V (1238.8, 1242. 8^) are seen shifted to lower 
wavelengths indicating stellar mass loss. The stellar 
material will be stopped in the H II region by coulomb 
forces, and it is suggested that most of the C+ ions in 
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I o INTRODUCTION 

One of the prime objectives of sounding- rocket spectroscopy 
has been to obtain ultraviolet interstellar line spectra which 
arise in the abundant species of C, N, 0, Si, S and Fe. Such 
observations require a hot background star which will provide 
sufficient flux against which the interstellar lines may be 
seen, and ideally, instrumental resolution of 0.l8 resolution 
or better. Perhaps the most satisfactory star for this purpose 
is C Oph. It is a bright, 09.5 V star with a large rotational 
velocity equal to 290 km sec -1 (Slettebak, 1956). The desirable 
consequence is that the confusion between stellar lines and 
interstellar lines is minimized. Furthermore, a recent and 
complete analysis of the visible interstellar spectrum of 
C Oph (Herbig, 1968) is available for comparison purposes. 
According to Herbig the cloud structure is comparatively simple 
consisting of a strong -15 km sec -1 component and a weak 
-29 km sec” 1 component. Under high resolution, however, these 
components show structure and additional components appear. 

The star is moderately reddened (Eg y = +0.32), and it exhibits 
a rich molecular spectrum in the visible spectral range. Thus, 
it is reasonable to expect that ultraviolet molecular band 
systems such as those of H 2 and CO will appear in the spectrum. 

An ultraviolet spectrum of C Oph from 1150 to 1700^ with 
the moderate resolution of 0.58 has been obtained by this 
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laboratory. Despite the fact that the planned for resolution 

of 0.1$ was not attained interstellar lines due to 0°, C + , 

*{* * 4 * 

0 5 Si a and S are clearly present as are those of the fourth 
positive system of CO which have been reported earlier (Smith 
and Stecher, 1971). The analysis of atomic and ionic lines 
constitutes the subject matter of this paper, and although 
the accuracy of the results is poor due to the limited 
resolution it is possible to reach some interesting tentative 
conclusions. 


H 



II. EXPERIMENT 

The instrument consists of two one-meter Wadsworth- type 
spectrographs in which spherical gratings serve as the objective 
as well as the dispersive components. One of these, intended 
to record the C Oph spectrum from 900 to 1450$, did not operate 
properly. The other provided an unwidened spectrum extending 
from 1162 to 1710$ with a nearly linear dispersion of 4.51$ mm 
and a resolution of about 0.5$. The spectrogram was recorded 
on Kodak 101-01 film which was exposed for 188 sec at a mean 
altitude of 144 km. Subsequently, the data were reduced to a 
relative flux distribution with wavelength, the wavelength 
scale being determined by a quadratic fit of observed features 
to the stellar lines of C III (1247.37$), Si IV (1402.77$), 
and He II (1640.43$). Through consistency checks made with 
identified lines the accuracy of the assigned wavelength scale 
should be better than ± O.IA between 1247 and 1640$. 
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III. DISCUSSION OF TOE DATA 


Segments of the reduced data are presented in Figure L 
In each segment a wavelength scale of l8 div” 1 is plotted on 
the axis of the abscissas, and a relative flux scale is plotted 
on the axis of the ordinates. The relative flux scales are 
not necessarily the same in each segment, and the axis of the 
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lines indicate the assumed value of the background stellar 
radiation. 

The first two segments show strong absorption lines with 

half-widths of aboht 0,68 which can be attributed to the 

triplet resonance lines of S II (1250.5, 1253.8, 1259,58). 

On close inspection the lines show a shift toward shorter 

wavelengths of approximately 0.l8 from the laboratory values 

-1 

which corresponds to a radial velocity of -34 km sec , This 
shift is particularly evident for the 1253,8 and 1259, 58 lines. 
Apart from a possible contribution to the 1253. s8 feature from 
an interstellar P I (1253.78) line there are no other reasonable 
identifications for these lines. As is shown in Section IV the 
column density of S + ions’ calculated on the assumption that all 
of the absorption arises either in one cloud or in two is far 
greater than would be expected on the basis of the measured 
column density of hydrogen atoms and an assumed solar S/H 
abundance ratio. It becomes plausible to assume, at least 
temporarily, that the S + ions are more widely distributed in 
.space and velocity than the double-cloud model would imply. -Such 



an assumption would help to explain the shortward shift of the 
observed lijies, and would make possible the derivation of 
considerably smaller S ion column densities on the basis of 
the measured resonance line equivalent widths. 

None of the other lines shown in Figure 1 are completely 
symmetric about their minima , but with the exception of the C I 
(1280. lS) and C I (1328. 8$) blends their minima show good 
correlation with laboratory wavelengths. It is difficult to 
say whether the difference between the implied radial velocities 
of the C° and S + particles on one hand and the C + , 0° and Si + 
particles on the other reflect a real difference in the location 
and velocities of the particles or is simply the result of 
differing strengths of the components which contribute to the 
observed feature. In this regard it will be shown in Section IV 
that the equivalent widths of the C II, 0 I, and Si II lines 
lie on the square-root part of the curve of growth, while those 
of C I and S II lie near or on the flat part of the curve of 
growth. Perhaps all that can be said is that there is probably 
some structure in the C I and S II lines related to the velocity 
and spatial distribution of the C° atoms and S + ions which 
will be revealed only when a spectrogram of greater resolution 
is obtained. 

The stellar spectrum itself must contain to some extent 
the resonance lines of the singly ionized species. Thus, the 
broad dip in the spectrum delineated by the dashed line and 
centered on 1335. 2^ is to be associated with the stellar 
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C XI (1334. 5, j. 335.72) lines. These lines do not appear in 
the appropriate stellar model spectra of Bradley and Morton 
(1969) , and are presumably weak. Therefore, because of the type 
of star and because of the broadening of stellar lines 
(~l.s2) due to the high rotation of the star it seems reasonable 
to place the stellar background as shown. In the case of the 
C I (1328.82) blend it was difficult to fit a smooth background 
line to the gross features of the spectrum; a straight horizontal 
line was used instead. There is reason to suspect, however, 
that the feature centered on 13292 is due in part to Si"*" 1 * ions 
in the stellar atmosphere through the predicted Si III (1328.82) 
line, and the original expectation that this absorption feature 
was due exclusively to carbon atoms is consequently in doubt. 

Likewise, the asymmetrical feature near 12802 which is 
primarily due to a C I blend is probably also formed in part 
by a stellar Si XII (1280.42) line with a possible contribution 
from a stellar Fe IV (1280.42) line. A measurement of an 
equivalent width was not attempted in this case. No identifications 
other than a blend of C I lines have been made for the feature 
centered on 1277.22. The feature is obviously asymmetric, and 
to estimate the equivalent width only the long wavelength side 
of the feature was used. The area measured is indicated by the 
cross-hatching. 

The same procedure has been followed in the cases of 
the Si II (1260.42), Si II (1304.42), and 0 X (1302.22) lines. 

Xn the first case, the long wavelength side of the profile is 


S' 



influenced by the (9-0) vibrational transition of the fourth 

12 16 * 
positive system of C 0 (Smith and Stecher, 1971) and a 

Si III (1301.28) line. In both the asymmetric and the symmetric 

features the areas used in the equivalent width determinations 

% 

are indicated by cross-hatching. 

The La line is observed in absorption on the same spectro- 
gram as the lines of the heavier elements and is effectively 
due to interstellar hydrogen atoms. A rectified flux 
distribution in the vicinity of 1215.78 which includes several 
features found in the wings of the La line is shown in Figure 
2. The equivalent width and column density of hydrogen atoms 
are found by fitting a computed contour shown by the dashed 
line to the observed rectified profile as discussed previously 
(Smith, -1970). 

The equivalent widths of the lines and blends discussed 
above are listed in the second column of Table 1. The quoted 
uncertainties reflect only a maximum uncertainty in the level 
of the stellar background flux. 



IV. COLUMN DENSITIES 


In general, the curve of growth procedures described by 

StrSmgren (1948) were followed in order to derive column 

densities. In the case at hand, all the data points lie on or 

between 'the flat portion and the square-root portion of the 

appropriate curve of growth, and this necessitates a computation 

of a curve of growth for each line or blend of closely spaced 

lines as in the case of the Cl multiplets. 

To illustrate the procedure we show in Figure 3 the 

computed curves of growth for the C II (1334. 5$) and C I (1277.28) 

lines. Plotted on the abscissa is log t where t is the 

o o 

optical depth in the center of the line and is given by the 
expression 


o 
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In this expression N is the number of absorbers cm 8 f is the 
oscillator strength and b. is equal to (X/c) (2kT/Am )^ 2 0 As 

A. O 

usual k is the Boltzman constant, T the temperature, and Am Q 

the mass of the atom or ion. In this case we set b s (2kT/Am ) 

o 

1.5 km sec which is the value determined by Herbig (1968) 
from fitting the Ca II H and K lines to Stromgren T s curve of 
growth. Plotted on the Ordinate is log (W/b^) where W is the 
equivalent width of an absorption line formed toy N absorbers cm" 
The optical depth in' a line, T(X), was calculated using the 
Voigt-f unction together with t For blends t(X) was simply the 
sum of the t(X) s s for each line. From the optical depths 
residual intensities were computed and a subsequent integration 


2 
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produced the desired equivalent widths. Thus, Doppler line 
broadening and radiation damping were included in the curve 
of growth computations. The effects produced by radiation 
damping are of primary importance in analyzing the present data. 
All oscillator strengths were taken from the publication of 
Wiese, Smith and Glennon (1966) and Wiese, Smith, and Miles 


(1969). 

Two curves for each line are shown to illustrate the 


difference between the assumptions of a single-cloud model and 


a double— cloud model. In the latter case 83% of the absorbers 
were in one cloud and 17% in the other, the proportions being 
estimated on the basis of the relative strengths of the Ca II 
lines in the -15 and -29 km sec -1 clouds observed by Herbig 
(1968). The data points reduced to log (W/b^) with b =1.5 km sec 
are plotted on the single-cloud curves. The boxes surrounding 
the data points indicate the magnitude of the uncertainty due 


to locating the stellar background flux (vertical dimension), 
and the effect of the uncertainty in the oscillator strengths 
(horizontal dimension). The column density of C + ions in the 

2p° ground state corresponding to the horizontal position 

^ -+0 5 17 -2 

of the data point as shown is (1.4 _q] 3 ^ 10 cm * The location 

of the data point well up on the square root portion of the 


curve of growth suggests that the asymptotic expression giving 
the equivalent width of a line in terms of the square root of 
the absorber column density could be used. The appropriate 
equation as given by Stone and Morton (1967) is 


n 



( 2 ) 
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Me & u 

p 

where the transition occurs from the lower level L to the upper 
level u with statistical weight g u which is broadened by 
spontaneous decays to lower levels 1 with statistical weights 
g l0 The sum is carried out over all permitted transitions with 

1 < u. As a check on the column density of C + ions in the 

2 o 

P l/2 s ‘ tate determined by the curve of growth method the appropriate 

quantities were inserted in equation (2) with the result that 
17 -2 

N =« 1.4 x 10 cm , i.e. both methods produce the same result 

well within the limits imposed by the uncertainties. This 

calculation shows that for b = 1.5 km sec - * equation (2) will 

yield sufficiently accurate answers, and further, that a 

distinction between single and double clouds is not meaningful 

for the C II (1334.5, 1335. 7$) lines in view of the uncertainties 

in the equivalent widths. Using Figure 3 one can show that when 

the uncertainties are ignored the column density of C + ions in 
2 o 

the P^/2 s ^ a ^ e f° r the double-cloud model is a factor of 1.3 
less than in the case of the single cloud model. Herbig (1968) 
has derived a value for b of 2.4 km sec" from a fit of the 
Na I (3302.4, t 303.0$) and Na I (5890.0, 5895.9$) doublets to 
Str’dnigren’s curve of growth. If this value of the turbulence 
parameter is used the C ion column density is reduced by a 
factor of 2.5 and equation (2) is no longer valid. 

The C II (1335,7$) line which originates in the ground 
state P 3/2 level and the resonance lines of Si II (1260.4$) 
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Si II (1304.48), Si II (1526.78) and 0 I (1302.28) were all 
analyzed in the same way as the C II (1334.58) line using curves 
of growth computed for each line and a turbulence parameter 
equal to 1.5 km sec - ' 1 '. The resulting column densities for all 
these lines are listed in the third and fourth columns of Table 1. 

In each case the reduction in column densities resulting from 
a change of interpretation on the basis of a single-cloud model 
to a double-cloud model is relatively small, a factor of 2 or 
less, and within the listed uncertainties. 

Such is not the case for the C I and S II resonance lines. 

It seems likely in view of the shape of the spectral features 

identified with C I lines that several lines of each multiplet 

are present in a blend which is either intrinsic or induced by 

the instrumental profile. For an initial estimate of the strength 

of each line involved in the C I (1277.28) and C I (1328.88) 

multiplets it was assumed that the lower levels were populated 

by collisions and that the population of a given level i relative 

to the ground level o is given by n^/ n Q * (g i /g Q ) exp (-X^kT) 

where the n ± is the number density of atoms in energy level i, 

and X- is the excitation energy of level i. The carbon atom 

3 

has three fine structure levels in the P ground state configuration 
with excitation energies of 0.0000, 0.0020 and 0.0054 eV for 
which J - 0, 1, and 2 respectively. If T « 100°K the fraction 
of atoms in the J » 2, 1, and 0 levels respectively is 0.44, 0.39, 
and 0.17. With this information curves of growth for the C I (1277.2 
and C I (1328.88) multiplets and for both the single and double- 
cloud models may be computed. Although in the case of a blend 


13 



t loses its usual significance it may be defined such that the 
o 

linear part of a curve of growth computed for a blend of C I 
lines coincides with the linear part of a curve of growth 
computed for a single line. This normalizing procedure has 
been followed in plotting the computed curves of growth for the 
C I (1277.28) blend shown in Figure 3. The results of this 
analysis for both the C I (1277. 2$) and C I (1328. 8$) blends and 
for both the single and double cloud models are given in columns 
3 and 4 of Table 1. If all the carbon atoms were in one cloud 
the uncertainties permit a reduction in the listed column density 
corresponding to the C I (1277.28) multiplet by a factor of 45. 

If the carbon atoms are distributed between two clouds according 
to the recipe given above the resulting column density is 40 times 
less than in the single-cloud. case and the uncertainties could 
increase this factor to 200. Thus the column density of the 
carbon atoms is very uncertain and quite dependent on the model 
used for analysis. 

The S II resonance lines are puzzling; the equivalent 
widths of these lines are quite large as are those of the 
C I resonance lines , but the strength of the S II lines cannot 
be explained either by the experimental uncertainties or by 
assuming that the S + ions, are distributed as prescribed by the 
double-cloud model. Analyzing the data using the curve of growth 
method described above the column densities listed in Table 1 
for both the single' and double-cloud models were derived; they 
are much too large if one expects the interstellar abundances 

* 

to be approximately those found in the sun c For example , the 



the equivalent width of the S II (1259,58) line yields a S + 

19 -2 

ion column density of 1,4 x 10 cm The column density of 

hydrogen atoms as previously determined (Smith and Stecher, 1971 ) 

from the La line is (4.6 ± 1.5) 10 cm . Thus, the approximate 

sulfur to hydrogen abundance ratio averaged over the distance 

from the earth to the nearest boundary of the C Oph Strdmgren 

-2 

sphere is 3,0 x 10 which is to be compared to the solar 
value of 2.0 x 10~ 5 (Goldberg, Muller, and Aller, 1960). As in 
the case of the carbon atoms the uncertainties in the column 
densities of the S + ions in the double-cloud model are large 
due to the fact that the log (W/b^) values lie close to the 
saturated part of the curve of growth. However, the uncertainties 
only permit a maximum reduction in the measured sulfur to hydrogen 
abundance ratio to a value which is still a factor of 10 larger 
than the solar abundance ratio. 

It is possible to apply the doublet ratio method to the 
1250.5 and 1253.88 lines of S II since both lines arise from 
the same level. Use of this method implies that the equivalent 
width of the S II (1259.58) line is in error since the line 

has the smallest measured equivalent width and yet the largest 

+ 17 -2 

f-value. The resulting S ion column density is 3.7 x 10 cm 

-4 

which gives an abundance ratio of 8.0 x 10 . The derived 

turbulence parameter, b, is 6.7 km sec - '*' which is considerably 
greater than the largest value for b (2.4 km sec”*') found by 
Herbig (1968). The ratio of the equivalent widths of these 
two lines, however, is sufficiently uncertain to permit the 
larger column densities and small turbulence parameter determined 



by the first method, and therefore little insight is gained by 
using the doublet ratio method. 

In view of these results it seems likely that the observed 
S II features are due to an unknown number of component lines 
which the instrument is unable to resolve, but which unfortunately 
has a most critical effect on the derivation of the S ion column 

-J ~ ~ ~ — -<s — : 

UCAiOJ. JL ^ ® 

While the determination of the oxygen atom column density 

is accomplished by the methods discussed above, the question 

arises as to how much of the oxygen is located in the earth’s 

atmosphere above the rocket. The simplest approach in answering 

this question is to integrate a model distribution of oxygen in 

the earth's exosphere outward from an appropriate mean altitude 

to obtain a telluric contribution to the total oxygen column 

density. For this purpose the summer-time model with an exospheric 

temperature of 1500°K contained in the ”U. S. Standard Atmosphere 

Supplements” (1966) was used. The column density of all the 

residual atomic oxygen obtained when the integration starts at 

the mean altitude of 144.5 km achieved during exposure is 
16 -2 

6.7 x 10 cm Stone and Morton (1967) estimate that 62% 

3 

of this oxygen is in the ground state from which the O I 

(1303.2^) arises. Consequently, the telluric contribution 

16 -2 

becomes 4.2 x 10 cm which is 2% of the column density 
determined on the basis of a single-cloud model. 



V. LINE-OF-SIGHT INTERCLOUD GAS DISTRIBUTION 


In order to analyze the data further it is necessary to 
obtain a model distribution of the inter-cloud gas in both the 
H I and H II regions. The procedure discussed by Herbig (1968) 
is followed closely here where one of his models with the following 
characteristics is used. The radius of the H II region surrounding 
C Oph is 15 pc within which the density of hydrogen ions and 
electrons is equal to 3 cm" . According to Herbig (1968) 

* 

these values are compatible with the data afforded by direct 

photography of the emission region surrounding C Oph and radio 

observations. Outside the ionization front a shell of hydrogen 

— 3 

with a density of 3 cm and an electron to hydrogen atom 

density ration, e, of 6 x 10“ 4 extends to 20pc. From 20 to 

170pc, i.e. to the sun, e remains the same and the hydrogen 

-3 

atom density is 0.07 cm . The electrons are presumed to 

arise from the photoionization of C, Mg, Si, and Fe which are 

present with solar abundances. The average hydrogen atom 

density along the extended 20 to 170pc path is reduced somewhat 

—3 

from Herbig’s value of 0.1 cm , reflecting a heavy reliance 
on rocket results (Jenkins, 1969) for determining typical 
intercloud densities. The temperature inside the H II region 
is 10 4 (°lO, and outside the H II region is 10 2 f*K). This model 
corresponds to one of the three which Herbig (1968) discusses 
with the omission of an interstellar cloud in which most of 
the visible spectrum arises. 

The column densities of atoms or ions in these three 
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extra-cloud regions may be found by performing the integrations 
designated in equations (13a, b) of Herbig (1968). When the 
case of carbon is considered these equations become 


N(C°) - R(C) t n 3 (H) dr, and 


N(C*) - R(C) | n 2 (H) dr 

* vs 


(3a) 

(3b) 


where 


D = e 2 n 2 (H) + h^ntfl) + h^g, 


(3c) 


Here, n(H) is the hydrogen atom density and R(C) is the solar 
C/H abundance ratio. The quantity h^ is defined as: 

(4) 


h i 5 —c r n(e> " r 

1 n(C ) a 


•&» Q 

where n(C ) and n(C ) are the number densities of singly ionized 

and neutral carbon atoms respectively, n(e) is the number 

density of electrons, T is the photoibniza tion probability per 

atom sec"\ and a is the recombination coefficient in units of 
3 -1 

cm sec . The quantity h ^ is related to the number density 
ratio of doubly to singly ionized carbon, n(C )/n(C ), in 
the same manner. The integrations in equations (3a, b) are to 
be carried put over the range of distances appropriate to the 
three different regions. 

The value of T may be calculated from Herbig*s equation 
reproduced here as: 


T ® 10“®h“^ J ° a \ u dX . 

o 


( 5 ) 



where h is Planck's constant, a is the photoionization cross- 
section, u is the energy density, and \ Q is the threshold for 
photoionization. The wavelengths are in Angstrom units. The 
energy density is comprised of two components, one due to the 
general galactic background, and the other due to the star. 

The energy density of the galactic background was determined 
from the results of Habing (1968). In order to estimate the 
stellar component of the energy density a line blanketed model 
of Bradley and Morton (1969) with T g =» 28640 (°K) and log g ** 4.0 
was used. The stellar radius was assumed to be 8R Q , and the 
integrations in equation (5) started at 227.88, the onset of 
the He II continuum. The photoionization cross-sections were 
calculated for all atoms and ions with the exception of Si and 

•Jo . _ *Jo 

Si from formulas given by Seaton (1958) . For Si and Si 
similar formulas of Chapman and Henry (1971) were used. 

The recombination coefficient for atomic oxygen (final 
4 o 

state) with T = 10 ( K) was kindly provided by Dr. R. Kirkpatrick. 
Seaton’s (1951) values for atomic carbon at 10^(°K) and 10^(°K) 
were used, and recombination coefficients at these two temperatures 
for singly ionized carbon were calculated on the basis of formulas 
given by Silk and Brown (1971). Hydrogenic values were used 
for all other recombination coefficients required. Table 2 
summarizes the results. Typical values of T fpr each atom or 
ion of interest at varying distances from the star are listed 
in the columns. The recombination coefficients are listed at 
the end of the columns for which they apply. 
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The photoionization coefficients for C , Si , O , and 0 

are set equal to zero outside the H II region since their 

ionization potentials are greater than that of hydrogen, and 

the energy density of extreme ultraviolet radiation (\<200$£), 

X-rays, and cosmic rays can for our present purpose be ignored. 

Using the model parameters and the values of T and a, for 

n'lixC h typica JL V a axe ^/x eoeii teu jl xx X a U X Z , xne 1IXXC ^1*3 X lODS 

in equations (3a, b) were performed. The resulting column 

densities are listed in Table 3 for each of the three regions 

through which the line-of-sight passes. This model distribution 

of atoms and ions will be used in the following sections to 

correct the measured column densities for contributions from 

the intercloud medium. In some cases the corrections are small 

and errors in the model are unimportant. In other cases, 

4. 

particularly that of C , the corrections are large, and 

knowledge of the correct distribution of the intercloud gas is 

+ 

important. Consequently, the column density of C ions inside 
a cloud determined from the present observation is quite 
uncertain. 
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VI. 


POPULATION OF THE 2 P° 3 / 2 LEVEL OF THE C + 


IONS 


The strength of the C II (1335. 7A) line is remarkable in 
that it represents absorptions from the excited, 0.0079 eV , 

2p 3/2 fine structure level of the ground state which, therefore, 
must be strongly populated. This is borne out by the results 
presented in Table 1. The radiative decay of this level has 
been recognized as an important cooling mechanism for inter- 
stellar clouds (e.g. Penston, 1970, and references appearing 
therein), but such a large population of the excited level was not 

altogether expected. Another observation relevant to the excitation 
-H + 

of the C ions is that a similar level in Si ions at 0.036 eV is 
not appreciably populated since no obvious lines arising from this 


excited level appear in the spectrum. If it is assumed that the 
excitation of the C + ions is achieved by collisions with gas 

3 

particles then by setting the mean energy per gas particle, ^ kT, 
equal to the fine structure level excitation a lower limit to the 
gas temperature may be estimated; this limit is 61 K. In addition 
to the fact that the C II (1334.5, 1335. 7^) lines are narrow, i.e. 
close to the limiting instrumental resolution of 0.5$, the apparent 
absence of Si + ions excited to the 0.036 eV , 2 j> 3/ 2 leve l seems to 
rule out the possibility that the C II lines are formed in the 
stellar atmosphere. 

The density of the gas may be estimated if it is assumed 
that the excitation is accomplished through collisions of C 
•ions with hydrogen atoms and electrons in the case of an H I 
region and with protons and electrons in the case of an H II 
region. This neglects collisional excitation due to Hg molecules 
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for which cross-sections are unknown to the author,, The 

observed population ratio of the 2 P°^ 2 to 2p °/ 2 levels in 
+ 

C ions is related to the excitation and de— excitation rates 
by the equation 


Q( 2 i 2 ^ 


Q(g } 0 + i) 


, _+ 

" ^l/2 } 


(7) 


1 3 

where Q(^, -) is the total collision excitation rate, Q(|, 1) 


is the total collision de-excitation rate, and w (~, i) is 

da da 


3 1, 

2 ! 2 ‘ 

the radiative decay rate. The quantities N 7 (cj^ 2 ) and N 7 (C*^) 
are the column densities of C + in the 2 P°^ 2 and 2p i /2 levels 
respectively within a cloud, i.e 6 the observed column densities 
corrected for contributions from the intercloud medium. 

Electron cross-sections were computed from equation 
of Bahcall and Wolf (1968) which is reproduced here as 


<0 v> 


(3 » 63 x 10~ 8 cm 3 sec" 1 ) 0 


~T/2 


4 


g l/2 


( 8 ) 


Here, < ov > e is the value of 0v averaged over a Maxwellian 
velocity distribution, |[) is the transition collision 

strength, T 4 is the temperature in units of 10 4 (°k) and S^ 2 
is the statistical weight of the 2 P ^y 2 ^ eve ^» For upward 
transitions the right side of equation (8) must be multiplied 
by exp(-X/kT) where X is the excitation energy 0 From the 
principal of detailed balance the downward rates may be obtained 



by dividing the right side of equation (8) by 53 / 2 ^ 1/2* For 
the quantity 0 (~, §) a value of 1.44 was used which was taken 

from Table 4 of Bahcall and Wolf (1968). 

Calculations of the excitation cross-section due to 
collisions with hydrogen atoms at various center-of-mass energies 
have been made by Weisheit and Lane (1971). It was necessary 
to average their computed values of ?v over a Maxwellian velocity 
distribution, and this computation was complicated by the necessity 
of obtaining the threshold value of o through extrapolation. 
Nevertheless, the value of <^v> for upward transitions was 
obtained in this way, and for downward transitions the upward 
rate was divided by (s 3/2 /s l/2 ) exp <-X/kT) . The results of 
these calculations are summarized in Table 4 which contains 
the collision rates per hydrogen atom due to both electrons 
and hydrogen atoms. In order to relate the electron collision 
rate to the hydrogen atom density, we have assumed that the 
ratio of the electron density to the hydrogen atom density is 

6 x 10" 4 . 

Now, to obtain the ratio of N' (C 3 / 2 ^ N ^ C l/2^ it is 
necessary to correct the derived column densities appearing in 
Table 1 for contributions from regions outside the cloud. A 
preliminary calculation of collision rates shows that there is 
no substantial excitation of the J =« 3/2 level outside the cloud 
in either the H II or H I regions since the densities in these 
regions are far too low. Thus, no correction need be made in 
the case of the excited ions, thus N^Cg^) ■ N ( C 3/2^* Making 



use of the column densities of Table 3 and the derived value 
for N(C^ /2 ) in Table 1 we find that N (C* /2 ) = 0.44 x 1G 17 C m‘ 2 . 
The ratio, N'CCg^^/N^CC*^) , is therefore equal to 2„3. This 
ratio is not possible since the maximum possible population ratio 
of the fine structure levels occurs in the thermal equilibrium 
limit and is equal to S^/2^^^/2 ~ 1-^ we a{ *opt column densities 

which result from using the maximum value of the stellar background 
level the ratio ., N'CC^/N '(C^/ 2 ) , becomes 1.5 reflecting 
principally on the decreasing importance of corrections. This 
fact illustrates how critical the estimate of the C + ion column 
density outside the cloud is - both in the case of C Qph and 
in other cases where the Stromgren sphere contains significant 


amounts of material. It should be noted, in this respect, that 
another way of reducing the N (cJ/^/N (C* /2 ) x^atio' is by changing 
the interstellar gas model so that the H I region with n(H) = 0.07 


cm 


■3 


extends all the way from the ionization front to the earth 
excepting the cloud itself. In this case the population ratio is 
1 0 4, again the importance of the corrections to the N(e^ ? ) column 
density is demonstrated. 


Since the excitation or de-excitation rates are simply 
the product of the hydi^ogen atom density and the collision 
rates per hydroghn atom equation (7) may be used to calculate 
n(H). Using a value for *(3/2, 1/2) of 2.36 x 10~ 6 sec" 1 
(Bahcall and Wolf, 1968) and the collision rates listed in 
Table 4 it is found for a population ratio of 1.5 that there 
is no hydrogen atom density sufficiently large to produce the 

oar 

required population of the 2 Pg/ 2 level at any of the temperatures 
included in' Table 4. If the population ratio is reduced somewhat 

2H 



arbitrarily to 1.0 then the derived hydrogen atom density 
becomes 3.3 x 10 4 cm -3 at 150°K, 1.2 x 10 4 cm -3 at 200°K, 
and 0.82 x 10 4 cm -3 at 250°K„ 

4 * 

While uncertainties remain in the corrections to the 
column density and in the collision cross-sections the hydrogen 
atom densities derived above are probably minimum values. As 

-3 

such, they are considerably larger than the density of 730 cm 
found by Herbig (1968) for the same interstellar gas model as 
is employed in this paper. There is an additional problem 
concerning the implied electron density within the cloud. Using 

the assumed ratio of the electron density to the hydrogen atom 

-4 4 

density of 6 x 10 and a hydrogen atom density of 10 , the 

-3 

implied electron density is 6 cm „ On the other hand, Herbig 

-3 

(1968) gives an upper limit on the electron density of 1.9 cm 

which is based on the measured line strengths of the Ca II 

(3933.7, 3968.52) doublet and the absence of the Ca I (4226.72) 

line. More recently, Bortolot, Shulman, and Thaddeus (1972) 

have established on the basis of the calcium lines that the 

-3 

electron density lies between 0.06 and 0.4 cm 

Corroborative evidence for cloud densities appreciably 
less than 10 4 cm” 3 comes from the CO bands of the fourth positive 
system found in the C Oph spectrum (Smith and Stecher, 1971). . 
Although it was impossible in this work to resolve line structure 
within a band, the band strengths were large enough to indicate 
that each band was comprised of several lines. In fact, the 
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results were consistent with the supposition that the first 
three levels of the ground state were populated by cosmic 
blackbody radiation. When a curve of growth was computed 
assuming that the molecules were in thermodynamic equilibrium 
with an ambient gas at 20°K a satisfactory fit of the data 
points to the curve of growth produced a turbulence parameter 


„ ^ 1 j. _ r\ Ary n - 
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Tne xowest value of the turbulence 

-1 


parameter determined by Herbig (1968) was 0.85 km sec and 

was based on the CH + molecular spectrum. Thus s it is not likely 

that the CO molecules are in thermodynamic equilibrium with a 

o 

gas of temperature equal to or greater than 20 K, but this 

condition would certainly exist if the cloud density was 

m 4 -3 
10 cm 

A possible alternative to collisional excitation of the 
2 o "I 3 

?3/2 level in C ions is to suppose that coherent scattering 

of 1334 . 58 quanta within the interstellar cloud populates the 
2 

upper D 3/2 - Leve ^ a rate sufficiently high that subsequent 
decays to the ^ 3/2 4eve ^ produce the observed ground state 
fine structure level population ratio. The expression giving 
the absorption rate per ion is simply B(v)F(v) where B(v) 
is the Einstein absorption coefficient and F(v) is the stellar 

Assuming the cloud 
to be 15 pc from the star, and using the flux measurements 
of Smith (1967) for C Oph at 1376§ the stellar flux entering 
the cloud is estimated to be 2.54 x 10” ergs cm - sec (cps) 
<» value which includes a correction made for interstellar 


-2 -1 -1 

flux in units of ergs cm sec (cps) 



rate of population of the “P” /o level in this way is 0.0503 x 10 


reddening. As an estimate of the effect of scattering, the flux 

2 

is increased in the cloud by a factor of t where t is the 

o o 

optical depth in the center of the C II (1334. 7^) line. The 
o 

factor t gives the order of the number of scatterings before 
o 

a \1334.7$ photon escapes assuming a 100% albedo from grains 

in the cloud and no redistribution of photon frequencies during 

2 

collisions. The optical depth in question is 5.0, is 25, 

— 17 —2 — 1 _ i 

and thus F(v) becomes 6.35 x 10 ergs cm sec (cps) . The 

-7 -1 

absorption rate per ion is then 2.97 x 10 sec . The 
probability that a C ion excited to the ^3/2 leve * will decay 
to the ^3/2 level is and consequently the estimated 

2^,0 — - ... A Ar AA - - 1 /\*~6 

3/2 

sec -1 . The de-population of the 2 Pg/ 2 level b y the reverse 

process is about 1/2 the population rate. When the population 

—6 —1 

rate is compared to the radiative decay rate of 2.36 x 10 sec 

(Bahcall and Wolf, 1968) it is apparent that this population 

mechanism cannot work. Furthermore, a more accurate treatment 

can only reduce the population rate. 

Thus, it becomes evident that the easiest way to resolve 

the discrepancies implicit in the assumption that most of the 

detected interstellar material is in one cloud is to suppose 

that a large" part of the C II spectrum arises outside the strong, 

-15 km sec -1 cloud. The fact that unacceptably high densities 

+ 

are required by the observations if the C ions are located in 
an H I region suggests that the extra C ions are located within 


27 



the Stromgren sphere. If the hydrogen of the "interior'* cloud 
is assumed to be completely ionized and at a temperature of 
10^°K then the density required to produce the relative ground 
state level populations in C + may be calculated. The electron 
collision cross-sections are given by equation (8) , and the 
proton collision cross-sections may be estimated from a graphical 


v: j/JL ooc 
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Wolf (1968). Assuming, as before, that the population ratio 
2 o 2 o •{• 

of the P 3/2 to P i/2 ^ eve ^ C is 1.5 then the density of 
electrons or protons necessary to produce this ratio is 


220 cm 


-3 


a not unreasonable value, 


2„o 


The question arises as to why the comparable ^ eve ^ 

+ 

at 0.036 eV in Si ions is not also populated to an observable 

degree. The answer is that the radiative decay rate between 

the fine structure levels is comparatively large, i.e, 2.13 x - 

sec (Bahcall and Wolf, 1968). This figure together with the 

proton collision cross-section, again estimated from Bahcall 

-3 

and Wolf (1968) , and the density of 220 cm derived from the 

—3 

C II lines combine to give a population ratio of 3.4 x 10 

The electron collision cross-sections are small compared to 

the proton cross-sections, and may be neglected. Thus, one 

2„o 


P 3 level in Si 


should not expect to see excitation of the 

— 3 

ions at densities on the order of 220 cm 

3 

For the case of the fine structure level in oxygen 

atoms an upper limit to the proton collision de-excitation 
cross-section is given by equation (31d) of Bahcall and Wolf (1968) 


zz- . 



reproduced here as 


<a(j = 1 to J = 2) v> =* (2.5 x 10 ' 8 cm 3 sec - 1 )<r 2 > (9) 

2 

where <r > is the mean square radius of the atom*s outer shell 

of electrons in units of 8 2 , and in this case takes on the 

value of 0.55. Multiplying equation (9) by g 1 /g Q exp ( -X/kT) *» 2.93 

to obtain the excitation cross-section, using a density of 

-3 3 3 

220 cm and a radiative decay rate between the P^ and P 2 

-5 -1 

fine structure levels of 8,95 x 10 sec (Bahcall and Wolf, 

1968) an upper limit to the population ratio of these fine 

structure levels is 0„10„ Thus, the longward side of the 

* # 

spectral feature centered on X 1304,4 a and due mainly to the 

3 

Si II (1304.48) line may reflect the population of the P^ 

level in o&ygen atoms located inside the StrUmgren sphere. 

In exactly the same way as outlined above it can be shown 

that the relative population of the P 2 level (0.0163 eV) 

3 + 

to the P Q ground level in N ions is 0.044. Thus one can 
expect to find an interstellar feature at X 10848 composed of 
a resonance line, N II (1084. 08 ), and two lines arising from 
the 3 P 0 level, N II (1085.7, 1085.58). The population of the 

3 

third fine structure level, P^, should be about a factor of 

3 

4.3 smaller 'than the P 2 level. Unfortunately, this region of 
the spectrum was not recorded due to a payload malfunction. 

Such an observation would provide useful evidence bearing on 
the hypothesis of collisional excitation of some interstellar 
atoms and ions. 

<>ir 

Assuming that the population ratio of the ground state 
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fine structure levels is 1.5 the column density of C ions 

17 -2 

inside the StrSmgren sphere is 1.67 x 10 cm . This result 

also depends on the assumption that all of the C II (1335. 7^) 

line arises in the StrSmgren sphere. The column density of 

15 

sodium in both atomic and singly ionized form is 1.33 x 10 era 
assuming solar abundances. If the material is located 4 pc 
from the star the photoionization and recombination coefficients 
for atomic sodium may be obtained from Herbig (1968) , page 255, 
Table 5 S and these values together with an electron density 
of 220 cm may be put in equation (4) to obtain the ratio of 
n(Na + )/n(Na°) which is equal to 132. Therefore, the column 
density of sodium atoms in the Stromgren sphere becomes 1.0 

__ I O Q 

x 10” cm which is approximately 1/5 of the sodium atom 

column density observed by Herbig (1968), Since in any case, 

the contribution to the observed Na° column density from regions 

outside concentrations of interstellar material is small the 

hypothesis put forward here that there is some material with 
' -3 

density equal to 220 cm within the Stromgren sphere is not 
basically inconsistent with the visible observations. 
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VII. MASS LOSS FROM C OPH 

It is possible that the hypothetical interior cloud is 
simply an ordinary H I interstellar cloud which was overtaken 
by the ionization front surrounding C Oph. Another possibility 
is that the interior cloud is formed by matter escaping from 
C Oph itself, and the evidence for this conjecture is now discussed. 

In Figure 4a, b are shown two portions of the C Oph 
spectrum containing evidence for what is interpreted to be stellar 
material moving radially outward from the star. Figure 4a shows 
two sti'ong features at 1233.2 and 1237.38 which are correlated 
v/ith the N V resonance lines at the laboratory wavelengths 

1238.8 and 1242. 8$ respectively. A relatively weak feature 

at 1240 . 5$ has been tentatively associated with Mg II (1239.9, 
1240.48) resonance transitions and A1 IV (1240.2, 1240.88) 
subordinate transitions. There is an A1 IV (1237. l8) line which 
lies close to the 1237.38 feature, but throughout the entire 
recorded spectral range there is no evidence that the A1 IV 
spectrum is strong. It is, therefore, unlikely that the 
A1 IV (1237. l8) line can account for such a strong feature. 

Figure 4b shows a similar situation prevailing for c’ H * + 
ions. Th^ minima of the large absorption feature located at 

1541.8 and 1544.58 are correlated with the laboratory measured 
resonance lines at 1548.2 and 1550.88 respectively. It appears 
that the shape of the spectrum contour between 1545 and 15568 
indicates that there is emission associated with this feature 
as well as absorption. The feature has been noted and reported 
previously (Smith and Stecher, 1971). As is shown in this 
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reference the (0,0) vibrational transition of the fourth 
positive system in the CO spectrum would be seen if there 
were sufficient background continuum radiation at 15448. However, 
based on the strengths of the bands appearing in the CO spectrum 
this feature would be weak. In addition, there are subordinate 
Cl IV lines at 1545.2, 1549.2 and '1551.38 which are also 
expected to be weak. 

Finally, in Figure 4c are shown the strong features 
associated with the Si IV (1393.7, 1402.78) resonance lines. 

The minima of these features occur very nearly at the laboratory 
wavelengths, and there is some small asymmetry to the features 
which exhibit more absorption on the shortward side of the minima. 
There are, however, no large shortward shifted absorption 
features associated with these lines. Another CO line 
corresponding to the (5,0) band in the fourth positive system 
measurably affects the short wavelength side of the Si IV 
(1393.78) li ne, but the strength of the molecular band is much 
less than the strength of the Si IV line. 

The standard interpretation of Figures 4a and 4b, and that 
which is adopted here, is that the C and N ions are 
accelerated by momentum transfer from the radiation field 
(Lucy and Solomon, 1967). The fact that the Si ions do not 
seem to share this radial motion can be explained if one assumes 
that at altitudes where the momentum transfer occurs most silicon 
is more highly ionized, i.e. in the form of Si s and the 
Si ions are at lower levels where the atmosphere has not 
achieved any observable outward velocity. With this in mind 



one might also expect to observe the resonance lines of 0 VI 
(1031.9, 1037.68) shifted to shorter wavelengths if the 
wavelength coverage could be extended down to the region of 
IP at 1025.78. Certainly, such an interpretation is extraordinary 
for a star whose classification is 09.5V, but it is nevertheless 
considered the simplest explanation of the observations. • . 

One can easily calculate the average radial velocities ' 

of the C and N ions corresponding to the minima of the 

shortward shifted absorption features; these are 1340 and 1230 km 

sec 1 respectively. Both of these values exceed the escape 

velocity from the star of 901 km sec 1 calculated under the 

assumptions that the mass of Q Oph is 17 M 0 and the radius is 

8R q . The ions, however, are slowed down in their outward flight 

by coulomb forces exerted on them principally by the protons 

in the H II region. At a temperature of IO^Pk) the protons have 

6 —1 

a mean velocity of about 10 cm sec which is at least two 
orders of magnitude less than the C ions. Following Spitzer 
(1967, 1968) a slowing down time, t g , can be defined so that 
the average change of velocity per collision is w/t g where w 
is the initial ion velocity. The expression relating t to the 
initial velocity is . 


*s * 
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(n£ ) 
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where m and m f are the masses of the ion and field particle, 
i.e. proton, respectively. The diffusion constant A D is 


33 



A 
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8n e n.Z Z„lnA 
f f 

2 ~ 
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--3 Ze 

where n. is the density of field particles or 3 cm , — 

X o 

Z e 

and f are the electrical charges of the ions and field particles 
c 

in emu units respectively, and In A is a "cut-off" factor which 
accounts for the fact that beyond a limited range the coulomb 


Application of equations (10) and (11) show that for the C 

and N ++++ ions the slowing down time is 1.89 x IG^ 1 sec and 

1.63 x 10^ sec respectively. If an average stopping distance a 

1, is defined by the relation T ~ (w/2)t then for the C " and 

s 

N ions 1 equals 4.1 and 3.2 pc respectively. The suggestion 
is made that the material leaving the star slows down and stops 
in the gas of the H II region at distances up to 4 pc where 
the ions assume a distribution among ion states determined by 
the ionization coefficient, F, and the recombination coefficient , 

a s characteristic of the radiation field and electron density 

— 3 

at these distances. If the electron density were 220 cm 
and the temperature 10^°K)the predominant ionic states would 

* 5 ” 

be C s N j 0 , P , and S . In this picture a low density annulus 
is built up at large distances from the star which probably 
shares the expansion velocity of the surrounding gas in the H II 
region. The extent of the absorbing region along the line of 
sight may be easily estimated if it is assumed that the population 


ratio of the ^£ 3/2 to ^ 1/2 • 1 - eve ^ s in ^ is 1*5, that all of 
the C in the F 3/2 s ' fcate is inside the StrSmgren sphere, that 
the proton density is 220 cm” , and that solar abundances apply. 





The answer is that the line-of-sight dimension of the absorbing 


region is 0 o 47 pc., 



VIII. ABUNDANCES 


As discussed in Section IV the derived column densities 

-j- 

of carbon atoms and S ions are quite uncertain due primarily 

to inadequate instrumental resolution,, Estimates of the column 
| 

density of C and Si ions and oxygen atoms can be made 3 however , 
and 0-0 tti p a re d to the column density of hydrogen determined from 
the La line« To do this it is necessary to subtract from the 
measured column densities the estimated contributions from the 
H II region. These are composed of two parts , one, an ionized 

_3 

gas with n(e) equal to 3 cm , and the other s an ionized gas 

-3 


with n(e) equal to 220 cm" 


The contribution of the first 


part may be obtained from Table 3. The contribution of the 

4- 

second part can be estimated if, as before, all the C ions 

2 Q 

in the p o /2 s ' tate ar ® assumed to be inside the Strbmgren 

2 o 

sphere and that the population ratio of the level to the 

level is 1„5„ Turning to Table 1 for the column density 

level it is found that the number of 


2p° 

P l/2 

4 - 2 o 

of C ions in the p g /2 

+ . ■ 2 o 

C ions in the Pg /2 contained by the interior cloud is 

17 -2 

0,93 x 10 cm if the largest value permitted by the uncertainty 

in locating the background continuum is used. This value is 

appropriate for use in conjunction with the 1,5 population ratio. 

To estimate the contribution from the interior cloud to 
. “f* 

the Si ion and oxygen atom column densities the total column 
.density of C ions, 1.93 x 10 cm , was used together with 
solar abundance ratios. In the case of oxygen atoms the column 
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density derived in this way was reduced by an additional factor 
of 0.1 to allow for the population of the P-^ level as discussed 
in Section V. No correction for a possible contributioh to the 
0 I (1302.28) line from telluric oxygen atoms was made since 
the standard atmosphere estimate discussed in Section IV was 
negligibly small, and the excitation of the O I (1304.98) line 
has been assumed due exclusively to collisions occurring inside 
the Stromgren sphere. 

The results are given in Table. 5, The identity of the 

atom or ion and the line from which the column densities 

are determined a*re listed in columns 1 and 2 respectively. The 

column densities corrected for contributions from the H II 

region appear in column 3. The ratio of the corrected column 

20 

density to the measured column density of hydrogen (4.6 x 10 
cm -2 ) is listed in column 4, and the solar abundance ratio is 
listed in column 5. Again, the uncertainties appearing in 
column 3 are due only to the uncertainty in locating the 
background stellar continuum level. The Si ion abundance 
ratio results from averaging the corrected column densities 

**\ 

determined from both the Si II (1304,48) and the Si II (1526,78 
lin^. The column density associated with the Si II (1260.48) line 
became negative when the correction was applied, and therefore, 
was not used. This absurdity means that the inaccuracies of 
both the measurement and the correction procedures are too 
large to enable one to derive a significant abundance ratio 
from the Si II (1260.48) line. Unfortunately, the same lack of 
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precision affects the abundance ratios determined from the 
other lines; it is* perhaps, surprising that the agreement with 
solar abundance ratios is as good as it is. All that can be 
said is that the present results give -no reason to believe that 
the interstellar abundances of carbon, oxygen, and silicon 
depart radically from solar abundances of these elements,, 

The question arises as to whether contributions from 
within the Stromgren sphere could account for the anomalously 
large equivalent' widths of the interstellar Cl and S II lines. 

A glance at Table 3 will show that in the case of the carbon 
atoms such a contribution would have to come from the postulated 

interior cloud since the column density through regions where 

^3 2.3 <->2 

n(e) = 3 cm is expected to be only 3.2 x 10 cm” „ Assuming 

* 

that in the interior cloud the column density of C ions is 

17 “2 —3 ~~ 4 - 

1 0 93 x 10 --cm , n(e) = 220 cm” , and that T/a = 2,5 x 10 

(see Table 2) it is easily seen by using equation (4) that 

the columnar density of carbon atoms is 1.7 x 10 a This 

column density corresponds to log t « 2 . and referring to 

O' 

Figure 3 it can be seen that the corresponding value of log 

(W^/b^) is almost included within the error estimates of the 

observation when the single-cloud curve of growth is used. If 

the carbon atoms were located in the . postulated interior cloud 

and in one or more clouds outside the H IX region there would 

be no difficulty at all in reconciling the observed equivalent 

5 «=>2 

widths of the carbon lines with reasonably small (« 2 x 10 cm” 



The situation for the S + ions is not so clear. Assuming 
that solar abundances prevail the total column density of S + 

16 

ions within the Stromgren sphere is estimated to be 1.0 x 10 
—2 * 

cm . This figure is at least a factor of 10 smaller than 
the observed S ion column density even when observational 
errors have been assumed to exist. Therefore, it is necessary 

•f 

to postulate that the absorption features attributed to S ions 

must to a large extent arise in regions outside the H II 

region which itself makes a relatively minor contribution.. 

Finally, it is possible to compare the corrected column 

densities of C + ions and oxygen atoms listed in Table 5 with 

15 -2 

the CO molecule column density of 1.8 x 10 cm previously 
measured by Smith and Stecher (1971). The result is that • 
along the line-of sight and outside the Stromgren sphere 
approximately 4% of the interstellar carbon and 0.1% of the 
interstellar 'oxygen is in the form of CO molecules. 



IX . SUMMARY 


Presented in this paper are observations of interstellar 
o + o + + 

lines arising in C , C , 0 , Si , and S found in the ultraviolet 
spectrum of C Oph. In addition, the stellar absorption lines 
of C IV (1548.2, 1550.88), NV (1238.8, 1242. 88), and Si IV 
(1393.8, 1402.88) are reproduced, the first two of which are 
shifted to shorter wavelengths by about 68. Determinations 
of interstellar column densities have been combined with a 
model for the distribution of intercloud material with the 
following results. 

■ 5 * o 

1. The abundance ratios obtained by dividing the C ,• 0 „ 

+ 

and Si column densities extending from the boundary of the 

Stromgren sphere to the sun by the hydrogen atom column density 

are consistent with the corresponding solar abundance ratios, 

the agreement being within a factor of 10 o The results imply 

that the C and S atoms and ions are distributed .among several 

clouds, and if there is no stratification of the species 

+ o + 

represented here the C , 0 , and Si atoms and ions are similarly 
distributed. Lines arising in the latter group lie on the 
square root portion of the curve of growth and are relatively 
insensitive . to the distribution of these particles. Such is 
not the case for the C° lines which are probably unresolved 
blends of closely spaced multiplet members which lie on or 
near the flat portion of the curve of growth. The column 
density of C atoms is in this case very sensitive to the 
distribution of atoms in the line-of-sight . 



2 . 


The 2 P°/ (0.0079 eV) level in C + ions is strongly populated. 

O / 


If the maximum column densities permitted by uncertainties in 


the background stellar flux are used the population ratio of 

2 „0 i 1 -1 ^ 

1/2 


the 


level to the “P~ /0 level is 1.5; It is suggested 


3/2 


-3 


that the excited C + ions are located inside the Stromgren 
sphere where proton densities equal to or greater than 220 cm 
would permit the ^ 3/2 ^- eve ^ ke PoP u ^ a i e ^ by collisions. 

3. The shortward shifted absorption lines of C IV (1548.2, 
1550.88) and N V (1238.8, 1242.88) indicate that matter is 
escaping from C Oph. However, the coulomb forces between the 
escaping matter and the residual protons inside the Strttmgren 
sphere will stop the stellar matter well short of the ionization 
front, indeed, probably within a distance less than 4 pc from 
the star. The suggestion is made that upon recombination in the 


H II region, part of the escaping C ions become C ions 
in the 2 P° / state, and thus contribute to the large population 

O/ £»' 

observed for this level. 
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TABLE 1 


EQUIVALENT WIDTHS AND DERIVED COLUMN DENSITIES OF 
LINE IN THE INTERSTELLAR SPECTRUM OF C OPH 


Line 

Identification 

W 

($) 

1 Cloud 
(cm -2 ) 


2 Cloud 
* (cm -2 ) 

C I (1277. 258) 

0.241 ± 0.094 

c.5 !1;|) 

10 17 

(3.6 

tZI. O. 

- 2 . 7 ' 

10^ 

C I (1328.83$) 

0.195 ± 0.064 

<®- 2 -t:P 

10 17 

(3.2 

+21.4. 
- 2.8' 

10 16 

C II (1334.53$) 

0.254 ± 0.037 

W-* t° 0 p 

10 17 

(1.2 

+0.4. 

-0.4' 

10 17 

C II (1335.68$) 

0.235 ± 0.039 

a.® 

10 17 


10 17 

0 I (1302.17$) 

0.262 ± 0.046 

«- 2 3:S> 

10 “ 

(1.7 

+1.0. 

-0.6' 

10 18 

Si II (1260.42$) 

0.281 ± 0.047 

<°- 87 !o'.24 )1016 

(o- 79 !S:l3> 

10 16 

Si II (1304.37$) 

0.130 ± 0.045 

_ +5.5. 
(6.5 _ 3 pl ) 

10 16 

(3.7 

+5.0. 
- 2.2 ' 

10 16 

Si II (1526.72$) 

0.195 ± 0.057 

(6.8 t^6> 

10 16 

(3.8 

+3.3. 

-1.8' 

10 16 

S II (1250.50$) 

0.156 ± 0.040 


10 19 

(2.6 

+3.7. 

-1.7' 

10 19 

S II (1253.79$) 

0.167 ± 0.048 

< 2 - 9 tip 

10 19 

(1.9 

+4.1. 

-1;9 ; 

10 19 

S II (1259.53$) 

0.118 ± 0.037 


10 19 

(6.6 

+16.8. 
- 5.7' 

10 18 
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PHOTOIONIZATION COEFFICIENTS, T, AND RECOMBINATION COEFFICIENTS 
a, FOR VARIOUS DISTANCES FROM C OPH 
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TABLE 4 


COLLISION RATES PER HYDROGEN ATOM FOR TRANSITIONS BETWEEN THE 
GROUND STATE FINE STRUCTURE. LEVELS OF C + 


Temperature 

(°K) 

. J 

-J* 

<o v > 

Electrons 

, 3 ~L 

(cm sec ) 

Hydrogen Atoms 


.1 

Ci 

3 

o 

6 

1.50 (-10) 

5.4K-10) 

100 


3‘ 

2 

1 

2 

1 . 86 (-10) 

6.74 (-10) 

ft 


1 

5 

3 

" 2 

1 . 66 (-10) 

7 o 01 (-10) 

150 


3 

2 

1 . 
2 

1 . 52 (-10) 

6 o 44 (-10) 


1 

2 

3 

" 2 

1.67 (-10) 

8.05(-10) 

200 

• 

3 

2 

1 

2 

l o 32(-10) 

6.35(-10) 

250 

1 

2 

3 

2 

1.64 (-10) 

8.60 (-10) 


3 

2 

1 

2 

1.18(-10) 

6 , 19 (-10) 
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TABLE 5 


ft 


ABUNDANCE RATIOS DETERMINED FROM THE INTERSTELLAR 

SPECTRUM OF C OPH 


Ion or 
Atom 

Absorption 
Line (8) . 

Corrected Column 
Density (cm -2 ) 

Abundance 

Ratio 

Solar 

Abundance 

Ratio 

c + 

1234.5 

17 

0.4 x 10 x ' 

0.9 x 10~ 4 

5.2 x 10‘ 4 

0° 

1302.2 

(1.9+°;g)10 18 

41 x 10" 4 

9.1 x 10" 4 

Si + 

1304.4 

(5 * 2 -3;i )l<)16 7 : 

12 x 10" 5 

3.2 x 10~ 5 

Si + 

1526.7 

< 8 -*3:$ )xo 7. 




H? 


.Cm** 


Bahcall, J. N. , and Wolf, R. 'A. 196S, Ap. J. , 152, 701. 

Sortolot', V. J. , Shulman, S., and Thaddeus, P. 1972, to be 
published in the Ap. J. 

Bradley, P. T. , and Morton, D. C. 1969, Ap. J. , 156 , 687. 

Chapman, ,R. D. , and Henry, R. W. J. 1972 to be published in the 
Ap. J. , 173 . 

Field, G. B. , and Hitchcock, J. L. 1966, Ap. J. , 146 , 1. 

Goldberg, L. , Muller, E. A., and Alley, L. H. 1960, Ap. J. 

Suppl 0 , 5, 1. 

Habing, H. J. 196S, -Bull. Astr. Inst. Netherlands,- 19, 421. 

Herbig, G. H. 1968, Z.f. Astrophys., 68, 243. 

Jenkins, E. B. 1969, in "Ultraviolet Stellar Spectra and Related 
Ground-Based Observations" (International Astronomical Union 
■ Symposium No. 36) eds. L. Houziaux and.H. E. Butler (Dordrecht: 

Do Reidel .Publishing Company) . 

Lucy, Lo B„, and Solomon, P. M.- 1967/ A. J.,. 72, 310. 

Penston, M. V. 1970, Ap. J. , 162 , 771. 

Seaton, M. J. 1951, M. N. , 111 , 368. 

. 1958, Rev. Mod. Phys., 30, 979. 

Silk, J. , and Brown, R. L. 1971, Ap. J. , 163 , 495. 

Slettebak, A.., 1956, Ap. J. , 124 , 173. 

Smith, M. 1970, Ap. J. , 160 , 595. 

.Smith, A. M. , and Stecher , . T. ; P. 1971, Ap..J., 164 , IA3-. 

Spitzer, L,* i 1967, "Phys'ics of Fully Ionized Gases" (New York: 

John Wiley & .Sons) . 

> 1968, ‘"Diffuse Matter' 1 in Spas©" (New York; John Wiley & Soa), 

< • i • 

■ • 0 

Stone, M. E. , and Morton, D. C. 1967, Ap. J., 149 , 29. 

.... . * . . . i ' — — 

Strom gran , 3. 1948,' Ap. J., 108 , 212. . 



o 


a 


U. S. Standard Atmosphere Supplements, 1966 (Washington: U.S. 
Government Printing Office). 

Weisheit, J. C., and Lane, N. F. 1971, Phys. Rev. A, 4, 171. 
Wiese, W. L. , Smith, M. W. , and Glennon, B. M. 1966, NSRDS-NBS4 , 
1 , 1 . 

Wiese., W. L. , Smith, M. W. , and Miies, B. M. 1969, NSRDS-NBS 22, 

2 , 1 . 


L tZ 



FIGURE CAPTIONS 


Segments of the C Oph spectrum containing various 
interstellar lines. Abscissae, wavelengths in 
Angstrom units. Ordinates, flux in relative units. 

Dashed lines indicate the assumed background stellar 
flux, cross-hatching the area used in computing 
equivalent widths. 

Rectified flux distribution of the La line. Abscissae, 

wavelength in Angstrom units. Ordinates, flux in 

o 

relative units. Features due to N , Si , and H^O 
which modify the La line contour are indicated. 

Dashed curve is a computed profile fitted to the core 
of the observed line. 

Curves of growth for the C I (1277. 2$) and C II (1334. 5$) 
lines computed assuming that interstellar matter is 
distributed in one cloud or two clouds. Abscissae, 
log of the optical depth in the center of the line. 
Ordinates, log of the line equivalent width divided by 
the Doppler width. • 

Segments of the C Oph spectrum containing the stellar 
lines of (a) NV (1238.8, 1242. 8$) , (b) C IV (1548.2, 
1550.82), and (c) Si IV (1393.8, 1402.82). Abscissae, 
wavelengths in Angstrom units. Ordinates, flux in 
relative units. Vertical lines labeled with numbers 
indicate the position of the possibly identified 
lines listed in the box. 
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